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I. MOTIVATION 



A few years ago a calculation of the sum of van der Waals interactions for a purely 
dielectric cylinder in the dilute-dielectric approximation led to a surprising null result 
(This was first found by an unpublished calculation by the first author of the present paper, 
and later independently confirmed by calculations by Milonni and by Barton This 
unexpected finding produced a quest to calculate the corresponding Casimir energy with 
the aim of verifying the predicted equality between both quantities, which was recently 
established in Ref. However, a physical understanding of why this Casimir energy 
should vanish (as also does that for a dilute dielectric-diamagnetic cylinder remains 
elusive. In addition to the physical interest of this subject, we recall that the authors of 
Ref. 0] commented at some length upon the comparative advantages and shortcomings of 
Green's function formalisms and mode summation methods for the evaluation of Casimir 
energies. A formal relation between these two approaches was given in appendix A of that 
paper. 

The first procedure is essentially a calculation of the expectation value of the stress-energy 
tensor expressed in terms of Green's functions and their transforms. This approach has 
proven to be remarkably fruitful and enlightening from the perspective of physical interpre- 
tation, as the quantities computed bear close relationships to easily identifiable observables 
expressed in terms of sources and fields. This was the method applied in Ref. [H to the 
object of our present study, i.e., the Casimir energy for a dielectric cylinder with the light 
velocity different on the inside and the outside. 

In contrast, the mode summation method rests on the concept of zero point energy and 
its basic expression as an infinite sum of eigenfrequencies. Even though this sum is a rather 
abstruse concept, from a prosaic viewpoint it has the appeal of simplicity. Since these 
eigenfrequencies stem from classical problems and are in many cases already known (e.g. 
through textbooks like Ref. 0), the only remaining task is to perform their summation. 
Once this question has been mathematically posed, the answer sought after is rather easily 
within one's grasp. 

In view of these considerations, the reflections in Ref. [lj itself, and the examples offered 
by many related works, it is sensible to say that the two procedures may complement each 
other. This is especially so, since these calculations still present a number of subtleties, 
and it is to be hoped that approaching the problems from disparate viewpoints may lead 
to improved physical insight. Specifically, the well-known divergence difficulties, which have 
recently received much attention 0], encourage us to examine the problem anew. With this 
idea in mind, we shall tackle some aspects of the problem treated in Ref. J3j by a particular 
variant of the mode summation method. 



II. MODE SUM 

According to Refs. 0, II, the eigenfrequencies to m ,p,k z of the Maxwell equations for an 
infinite material cylinder of radius a, oriented along the z-axis, with permittivity and per- 
meability (eijfMi), surrounded by a medium with permittivity and permeability (£2,^2), are 
the zeros of the following function: 

f m {k z ,u) = 0, m = 0,±l,±2,..., k z eR, (2.1a) 
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f m (k g} u) = ^ 



a 4 u 2 k 2 

A™(x, y) A™(x, y) - m 2 ^^( £l/Ul - £ 2M2 ) 2 J^(x) ^(y) 



(2.1b) 



with the following abbreviations 



2i 

A = , x = Aid, y = A 2 a, 

IT f 

A m (x, y) = viy J m ( x ) H m {y) - /i 2 xJ m (x) H m (y), (2.2) 
&™(x,y) = £iyJL( x ) H m(y)-£2xJm(x)H' m (y), 

A 2 = EifiiU 2 - k z , i = l,2. 

As usual, m is the azimuthal quantum number, k z denotes the momentum along the axis 
of the cylinder, and p labels the zeroes of f m (k z ,u>). Here, for y > 0, H m (y) = Hm(y). 
Note that f m = — A~ 2 2, where S is the same denominator introduced in Ref. |J] and the 
A~ 2 factor has been introduced for convenience. Recalling that the velocities of light in 
each media are given by q = (e^/ii)™ 1 / 2 , i = 1,2, we see that Eqs. ()2.1aj) . (J2.1b|) exhibit 
a peculiar feature of this situation: if c\ ^ c 2 , the second term in (|2.1b|) comes into play, 
and it is no longer possible to decompose the mode set into zeros of A™ and zeros of A™; 
that is, the transverse electric and magnetic modes become entangled. When medium 1 is 
purely dielectric and medium 2 is vacuum, e\ = s, fi\ = 1, £2 = ^2 = 1, one may write 
uj = a^ 1 (y 2 + k 2 ) 1 ^ 2 and x 2 = y 2 + (e — l)(y 2 + k 2 ), where k = k z a. We will also set c 2 = c = 1. 
Formally speaking, the Casimir energy per unit length is given by the mode sum 



2 J-00 2tt 

m,p 



This quantity is related to the Casimir radial pressure Pc through 



Pc = —Sc- (2.4) 



na 2 



The sum in f!2.3|) is divergent and will be regularized by changing thepower of uj, which 
power becomes the argument of the associated zeta function. By Ref. |7( we already know 
that, up through the order of (e — l) 2 , there should be no ambiguities, i.e., no logarithms 
of arbitrary scales, because the heat kernel coefficient which would multiply them is of 
0((e — l) 3 ). Therefore, we are permitted to define 



v m,p ~ — ' m,p 



as a function of the complex variable s, without any additional mass scale. It is proposed 
to make sense of Eq. (J2.3|) by means of analytic continuation of this function to s = — 1, 
namely, 

E G = lim £ c (s). (2.6) 

s^ — 1 

For any given values of k, m, the remaining sum over p may be rewritten, with the help of 
the residue theorem, as a contour integral in complex y plane. Thus, 

fe w - r- 1 fli±^J c «yy fe 2 + ^ m /„ (2.7, 

J -°° m=-oo JL ' 
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where C is a circuit which encloses all the y values corresponding to the positive zeroes of 
f m . This approach is often referred to as the argument principle Quite often in the 
application of this method, an asymptotic form f mj3S of f m is found and then ln/ m as is 
subtracted from ln/ m in the integrand. Actually, the factors introduced in (j2.1b|) relative 
to the original f m of Ref. p| amount to dividing that function by the leading term of its 
asymptotic behaviour. At any rate, the full form of this asymptotic behaviour is presumably 
related to the limit of unbounded space, which is already available as the bulk contribution 
calculated in Ref. 0. 

Next, the dilute-dielectric approximation will be made by expanding the logarithm func- 
tion of Eq. ()2.7|) in powers of (e — 1), choosing y as the independent variable and taking x 
as a function of y and k: 



+ 



L° ml ( y ) + LUy)(y 2 + k 2 ) 
L™M + LZ(y)(y 2 + k 2 ) + LZ(y)(y 2 + k 2 ) 2 + L% 2 ( y )( y 2 + k 2 ) k 2 (s - i) 2 



e-l) 

2\ i r 20 



+ 0((e 



where 



(2.8) 



L lni(y) 
L™ 2 (y) 



T 10 

L m2 



(y) 



L 2 l(y) 



Ay 



yJL(y)H m (y), 



^°\y), 



y 2 Jm(v) H 'm(y)i 



2A 2 
1 



A 



2A 2 

L 2 ^(y) + L™(y), 



A^\y)J' m {y)H m {y) + - ( J' m {y) + y ( 1 - ^ ) Uv) ) H m {y) 



L 2 l A (y) 



L 2 l B (y) 



AAy 2 



AA 2 y 



Ag-°>(v) - -A^(y) 
y 

r(A^(y)) 2 , 



m 



A 2 y 4 

Here, the notation is 



Jm{y)H 2 m{y)- 



^\y) 



^'A m (x,y) 



for j 



1,2, 



(2.9) 



(2.10) 



SI =£2=1,^1 =fJ,2=l,X=y 



where A m stands for either A™ or A™ (in the free space limit there is no difference). 
Moreover, note that 



A m (x,y)\ £l -_ 



S2=l,fll=fl2=l,X=y 



-yW[J m (y),H m (y)] = A, 



(2.11) 



where W denotes the Wronskian. Explicitly, 



,(i,o) 



(y) 



i 



[y 2 J' m {y)H' m {y) + (y 2 - m 2 )J m {y)H m {y)] - {J m {y)H m {y))\ 



A%°\y) = (A£.°)(y))' 



m 2 + 1 

y 2 



A. 



(2.12) 
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where (A^'°' ) (y)Y = —A^,'°\y) ^ A^°\y), as shown. Once this task has been accom- 

dy 

plished, we insert Eq. ()2.8|) into Eq. ()2.7|) and perform the k integration. 

Consideration of formulas ()2.7|) - (j2.12|) shows the need of dealing with integrals of the 
type 

/OO f 
dk / dyyF{y){y 2 + k 2 )- a k 2 ^ (2.13) 
-oo JC 

where C is the contour specified above, and F satisfies F(—iv) = F(iv) for v G R, as well 
as having adequate asymptotic properties (see below). Looking at the (y 2 + k 2 ) powers in 
Eqs. (12 .7|) . (|2.8|) . we see that in the required cases a = s/2 + 1, s/2, s/2 — 1, and (3 — except 
for one integral with (5 = 1. Since a is just a translation of s/2, analytic continuation in s 
amounts to analytic continuation in a. 
Straightforward integration for k yields 

I = B (p + ±a-l3-^ J dyy 2 - 2 ^ F(y), (2.14) 

r(x)r(y) 

where B stands for the Euler beta function ~B(x,y) = — r-. Now we rotate the C 

T(x + y) 

contour so that it consists of a straight line parallel to and just to the right of the imaginary 
axis, closed by a semicircle of infinitely large radius on the right. 1 The branch line of the 
^2-2q+2/3 f unc tion, which starts at the origin, is placed so that the circuit does not cross 
it, that is, it lies along the negative real axis. (In the limit where the vertical part of C 
overlaps the axis, the origin could be avoided by a small semicircle, and, eventually, the 
integration along this infinitesimal part would vanish as s — > —1). Further, it is assumed 
that, combining the asymptotic behaviour of F and the possibility of varying the a value 
as necessary, the contribution from the large semicircle vanishes when its radius tends to 
infinity. Thus, the y integral reduces to an integration along the vertical parts of C, where 
y = ±iv, and one has ^ 2_2a + 2 / 3 = e ±iK{i-a+f3) v 2-2a+2p on upper and lower segments. As 
a result, 

7 = -2iB a -13- sin (ir (a-p- J dv v 2 ~ 2a+2/3 F{iv). (2.15) 

71 

Applying the reflection formula T(z)T(l — z) = to the gamma functions in the Euler 

sin7T2; 

beta function, we re-express I as 

1 = -2iB ((3 + ^,l-a\ sin(vra)^ dvv 2 - 2a+w F(iv). (2.16) 

Note that for the a values corresponding to s = — 1 (as = 1/2, —1/2, —3/2), and for f3 = 0, 1, 
the pole of the B function and the zero of the sine function in Eq. (j2.15j) combine to give a 
finite product, while in Eq. (J2.16|) each of these two factors is already finite. 



1 Another way to describe this rotation is in terms of a purely mathematical transformation, based on the 
required analyticity of the underlying Green's function. See Refs. 

HQ- 
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Applying formula (|2.16|) to Eqs. (|2.7jl . ()2.8|) yields the following result: 
£ c (s) = £ci(e - 1) + £ C 2{s)(e - l) 2 + 0((e - l) 3 ), 



where 



h sa s 1 / T s 
B 



2 2tt 2 ~U'~^) Sin HD E ./ 

ft sa s_1 
2 2vr 2 



m=— oo 



Tfl— — OO 



and 



/ 1 \ /*oo 





n 


sa s 1 


B| 




~2 


2tt 2 




4° 2 ( S ) = 


h 




B| 




~2 


2tt 2 


fS^W = 


ft 




B| 




~2 


2tt 2 




n 




B| 




~2 


2tt 2 



2J~- V2, — Jo 

' m=—oo 

)oo 
sin (-^) E 
Tfl— — OO 
OO ,. x 



j 4-s r 20AB/ • \ 



Tfl— — OO 



A. First order in (e — 1) 

Taking £ci( s ) from (|2.18|) . and L^ ll (zu ) from (|2.9|) one is led to 



s-l 



1 s 
2'~2 



(2.17) 



(2.18) 



(2.19) 



sUs) = ~ ^M^-iHH^) £ /_ ^^»^)- (2.20) 



The beta and sine functions are finite at s — —1. As for the integral, it will be rewritten by 
introducing the factor 1 = — vW[I m (v), K m (v)] = —v[I m (v)K' m (v) — I' m {v)K m {v)} for each 
m: 



POO v 

JO „ 



Tfl— — OO 



POO 00 POO 00 

- / dvv 2 ~ s J2 Uv)I' m (v)K m (v)K' m (v) + I dvv 2 ~ s J2 Ci^K 

m — — i-vi ~m — — rv~< 



(2.21) 



[v . 



We carry out the summation over m by using the addition theorem for the modified 
Bessel functions: 



J2 I m (kr)K m (kp)e m * = K (kR(r, p, </>)) 

ra=— oo 

i?(r, p, 0) = a/V 2 + p 2 — 2rpcos0, p > r. 



(2.22) 
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Expressions for products of four Bessel functions are found by multiplying differentiated 
versions of the Bessel function addition theorem (|2.22|) . integrating over 0, and setting 
kr = kp = v. A change of variable u = sin ~ is made. Next, instead of proceeding with the 
evaluation of the resulting u integral, the expression is multiplied by a v power involving s, 
and integrated over v, with the help of 

Y ^ i~x+tj,+u ^ y ( i—\-i-i+v ^ y ( i—x+fi— y ( i—x—fi—u >^ 



dxx~ x K^(x)K u (x) = 2 



-2-A J 



^2.23) 



r(i-A) 

(from result (6.576.4) in Ref. HI). Then the u integration is carried out afterwards. Thus, 
the emerging expressions depend on s only through gamma functions (see also Refs. [H, 0]). 
In fact, the results in formula? (81) of Ref. Q may be viewed as intermediate steps, the final 
result being for example 



dvv 



dvv 



2-s 



m=-oc 

<x> 

E I m {v)I' m {v)K m {v)K' m {v) 



m=— oo 

i r(^)r»(^)r(^) r(§) 



m=— oo 



r(3-a) 



r(f)' 

(2.24) 



Even if the left hand side of each such integral is not initially defined for s = —1, the 
right hand side provides the desired extension to this value through the existing analytic 
continuations of the gamma functions themselves. Moreover, every result displays a single 
pole at s = — 1 in the gamma function in the last factor of each denominator, while the rest 
of the expression is finite at this point. Therefore, each expression has a zero of order one 
at s = — 1. 

Because Eqs. (I2.24j) show that the two integrals in (|2.21|) have the same value (even before 
setting s = — 1), we conclude, in a neighborhood of s = — 1, 



£%As) 



0. 



-ci\°) — w - (2.25) 

Equation (J2.18)) shows that the £(ji( s ) contribution involves the integration of the L l ml {iv ) 
function, determined by Eqs. ()2.9|) . (j2.12|) . In consequence, 



i . . h sa s 1 „ 

eL( 8 ) = =- b 

C1K J 2 2vr 2 



m=— oo 



dvv 



2-s 



2 2 

2 



sin 7T 



m 



I' m (v)K' m (v) -[! + —) I m (v)K m (v) + ~{I m {v)K m {v))' 



[2.26) 



Bearing in mind the useful property 
d 



dv 



v 2 r m (v)K' m (v) - (v 2 + m 2 )I m (v)K m (v)] = -2vI m (v)K m (v), 



(2.27) 



we apply partial integration omitting a 'boundary term' which vanishes for a given s range 
(2 < < 3). Doing so, we find that Eq. ()2.26|) becomes 

,8-1 



t ^ S) ~ 2 2tt2 



B 



2' 2 



sin 7T 



x 



/ dw 1 - 8 E (im(«)^m(«)y+i / dvV 2 ~ S E 

Jo 1- ~ s Jo 

m= — OQ nm — — rv-* 



m=—oo 



[2.28) 
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Because of this s restriction, the integrals in ()2.28|) cannot be directly taken at s = — 1. 
However, if this difficulty is disregarded, we may formally set s — — 1 and find 



= -1) -- — 2 / dvy2 Yl (Im(v)K m (v))' - - — 2 / dvv 3 S~] I m (v)K m (v). 

v m=— oo v m=— 00 

(2.29) 

One could argue that the first part can be dismissed as a mere contact term (since it may be 
shown from Eq. ()2.22|) that it is local in v). The second part of ()2.29|) gives a contribution 
which exactly cancels the bulk contribution found in Ref . , where the same type of formal 
expressions was employed. [See formulas (72), (78) there. In that paper what was evaluated 
was the Casimir radial pressure Pc, which is related to Sc through Eq. ()2.4j) .] 

Perhaps a better argument is to again use the technique of multiplying each term in the 
m summation by 1 = — vW[I m (v), K m (v)], and turn the initial expression into a linear com- 
bination of integrals with summations of products of four Bessel functions [like Eq. (12.24)) ]. 
That linear combination would yield an identically null result — one that is zero for any s 
value — by virtue of the symmetries seen in Eqs. (|2.24|) under interchange of different Bessel 
function types (see also comment after Eqs. (80) in Ref. j4|). Hence, near s = — 1, 

S l cx {s) = 0. (2.30) 

From another viewpoint, according to the reasonings in Ref. ^3| (and references therein), 
which dealt with a similar problem for a dielectric ball, all linear terms in (E2 — £1) have 
to be subtracted because they represent the self-energy of the electromagnetic field due to 
polarizable particles. Therefore, one simply must remove the linear part, no matter what 
precise form it has. This, of course, is the physical basis for removing the bulk energy 
contribution. 



B. Second order in (e — 1) 

We start with the part called £c^(s), because its evaluation is most similar to that of 
the Sq^s), £ci(s) contributions. Using the fourth line in Eq. (j2.19|) . the fifth in Eq. (|2.9|) . 
expressions f!2.12|) with y = iv, introducing, again, a unit factor in terms of the Wronskian, 
and following the same argument that led to (|2.30p . one arrives at 

£%*(8) = 0, (2.31) 

near s = —1. Next, selecting the lines in Eqs. ()2.19|) . which define £^ 2 (s), ££? 2 (s), S^i 8 ), 
S&(s), the parts of Eq. (El which determine L™ 2 {y), L^° 2 (y), L™g{y), L^ 2 {y), the form of 
Am°'(!/) given by (j2.12|) (its square for the case of L^ 2 B (y)), and taking imaginary arguments 
y = iv, one finds 

S °c^ = l^ B ^-^ sin (-4) r^ 2 '* £ C{v)Kl{v\ (2.32a) 



S 



£l\{s) 



h sa s 1 ( 1 s\ . ( s 

— B -, 1 sin 7r- 

2 4vr 2 \2' 2 \ 2 



x / dvv 
'o 



2-s 



E 



m 



-(« + —] lL{v)K m {v)K' m {v 



(2.32b) 



£ 2 f( S ) 



ft 












2 




v2 




POO 






x dv 


V 2 




Jo 





2-s 



S 

~ 2 

oo 

E 



Tl 2, 



«)JC(«) + Il{v)K' m \v) 



m=— oo L - 
.2 ™2 



2(1 _ _ m 2 )I^)/>)K m (<;)<>) 



+ v 2 I'J{v)K' m \v) + [v 2 + 2m 2 + ^\ P m {v)Kl{v 



m 



+ 2vI' m 2 {v)K m {v)K' m {v) + 2v I m {v)I' m {v)K' m \v) 



■I 2, 



m 2 \ r , , „, , s „ , s „ / m 2 \ 2 



-I - V ) - 2 Iv + — ) I 2 m (v)K m (v)K' m (v) 



(2.32c) 



£&(s) 



h sa s 1 / 3 1 s\ . / s 

t- B -, l sin 7i- 

2 2tx 2 V2 2 V 2 



cfoi;- £ m 2 I 2 m (v)K 2 m (v). 



m=—oo 



(2.32d) 



The f integrals are again evaluated in the manner illustrated in Eq. (|2.24|) . All of the 
resulting formulas exhibit zeros of order one at s = — 1, —3, —5, .... Since all the beta and 
sine functions in Eq. (I2.32|) are finite at s = — 1, each of these quantities will just yield 
(D(s + 1). In fact, after calculating the coefficients, 

' (s + l) + C((s + l) 2 ), (2.33a) 



'C2\ 

£™(s) 



h 



192 7T a 2 

h 2sLr^ + v + °« s+1 ^ 

h _„„ l J s + l) + 0((s + l) 2 ). 



h 



(2.33b) 
(2.33c) 

■■('•2\-"i - "Tp^r, ^7; ■ ; v, i .i ; /. (2.33d) 

2oU4 Ti or 

From Eqs. ()2.3H) and (12.331) . we see that lim £02(3) = 0, i.e., the (e — l) 2 contribution to 
the Casimir energy per unit length in the dilute-dielectric approximation is zero. 

III. CONCLUDING REMARKS 



Applying an analytic regularization which changes the eigenmode power, we have calcu- 
lated — to the order of (e — l) 2 — what in Ref. [3J] is called a pure Casimir term, i.e., the 
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convergent component of the energy depending only on h, c, the electrostatic polarizability of 
the material and the dimensions of the body. Using the words of Ref. a forest of gamma 
functions has grown out of an analytic continuation. Although the applied technique might 
be regarded as somewhat physically opaque, its relations to more transparent regularizations 
have already been studied (see e.g. Refs. [HI). I n Ref- 0] the only relevant contribution 
of the bulk part took place at the order of (e — l) 1 . Since in the present regularization 
the corresponding term vanishes by itself, to omit or include a separate bulk part leaves 
the outcome unchanged. This remark is true as far as the analytic regularization method 
proposed in Sec. 6.1 of Ref. [H is concerned, but not so for the (regulated) numerical method 
presented in Sees. 6.2 and 6.3 there, where a detailed numerical cancellation occurs between 
the bulk and cylinder parts. As a result, the pure Casimir term is seen to vanish to order 
(e — l) 2 . Comparing this calculation to the one offered in Ref. a decrease in difficulty 
can be appreciated in the derivation of the expression to be evaluated. However, we should 
stress that a substantial part of this merit is not in the nature of the method itself, but in 
the use of a previously known equation for the classical modes of this problem. Another key 
ingredient has been the exploitation of a Bessel function addition theorem 0, 0, 0] . 

Ref. exhibits the presence of a divergence at the order of (e— l) 3 . Among all the possible 
divergences which show up by point-splitting or ultraviolet cutoffs, the ones which survive 
in analytic regularization schemes introduce an unavoidable ambiguity parametrized by an 
arbitrary length or mass scale. This issue is viewed with more or less concern, depending on 
the authors (see comments in Ref. 0). The standpoint of Ref. was to admit that the 
problem is simply ill-defined because its posing constitutes an idealization (the permittivity 
e as a function of the radial coordinate is treated as a step function). 

In any case, the reason why the Casimir energy term of order (e — l) 2 is zero for the 
cylinder, while it is finite for other geometries, remains rather mysterious. It is clear that 
we still have some way to go to understand quantum vacuum energies. 
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